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SOME ALGEBRAIC IDENTITIES FOR THE a-PERMANENT 



HARRY CRANE 



Abstract. We show that the permanent of a matrix is a linear combination of determinants of 
block diagonal matrices which are simple functions of the original matrix. To prove this, we 
first show a more general identity involving a-permanents: for arbitrary complex numbers a 
and fi, we show that the a-permanent of any matrix can be expressed as a linear combination 
of jS-permanents of related matrices. Some other identities for the a-permanent of sums 
and products of matrices are shown, as well as a relationship between the a-permanent and 
general immanants. We conclude with a discussion of the computational complexity of the 
a-permanent and provide some numerical illustrations. 



1. Introduction 
The permanent of an n X n C-valued matrix M is defined by 



(1) perM ^YjH 



M 



jMi)' 



aeS„ /=! 



^ ■ where Sn denotes the symmetric group acting on [n] := {1,. ..,n}. Study of the perma- 

nent dates to Binet and Cauchy in the early 1800s [11]; and much of the early interest in 
[-^ ■ permanents was in understanding how its resemblance of the determinant, 

^ : (2) detM := Y^ ^S^i'^) [1 ^>'^(/> 

where sgn(a) is the parity of a € 5„, reconciled with some stark differences between the two. 

An early consideration, answered in the negative by Marcus and Mine |8|, was whether 
k^ ■ there exists a linear transformation T such that per TM = detM for any matrix M. In his 

j_j I seminal paper on the #P-complete complexity class. Valiant remarked about the perplexing 

. ^ .■ relationship between the permanent and determinant. 

We do not know of any pair of functions, other than the permanent and 
determinant, for which the explicit algebraic expressions are so similar, and 
yet the computational complexities are apparently so different. (Valiant |fT7[ , 
p. 189) 

In this paper, we hope to give some insight to Valiant's remark, by bringing forth the 
simple identity 



(3) (-1)" perM = Y, (-1)^*" det(M • n), 



nefi 



["] 
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2 HARRY CRANE 

which expresses the permanent as a linear combination of determinants of block diagonal 
matrices. In (HJ), the sum is over the collection ^[„] of set partitions of [n] :- {1, ... ,n},#n 
denotes the number of blocks of ti € 'Pin], x^^ :- x{x - !)■ ■ -{x - j + 1) =: (-x)^-' (-!)•', and 
det(M • n) :- Ylben detM[&], a product of determinants of the submatrices M[b] with rows 
and columns labeled by the elements of each block of n. Equation (O is an immediate 
corollary of our main identity (O for the a-permanent II18| w^hich, for any a € C, is defined 
by 

I! 

(4) per„M:=2]^a*'^[]My,,o), 

aeS„ 7=1 

w^here #o denotes the number of cycles of a e S„. The a-permanent generalizes both 
the permanent and the determinant: perM = per^M and detM = per_j(-M); and, when 
M = }n, the nxn matrix of all ones, IHl coincides with the generating function of the Stirling 
numbers of the second kind: a^" = L^'=o s{n,k)a^'^, where s{n,k) := #{n e P[„] : #ti - k]. Our 
main identity, 

(5) per^^ M^ Yj '^^^^ P^^ai^ ■ ") for ^11 a,peC, 

expresses the a-permanent as a linear combination of |S-permanents, for any choice of a 
and jS. This identity, and its corollaries, could be insightful to understanding the apparent 
gap between the computational complexity of (HJ and ^. We discuss these observations 
further in section |3.3[ and state a conjecture about the computational complexity of the 
a-permanent. 

In addition to (O, we show other identities for the a-permanent of sums and products 
of matrices. We separate these identities into two main theorems, calling the first the 
Permanent Decomposition Theorem. 

Theorem 1.1 (Permanent Decomposition Theorem). For any a, jS 6 C and M e C^", 



nePt 



["] 



where 'P[„] is the collection of set partitions of [n] :- {!,..., n}, p^J :- |S(j3 - 1) • • • (jS - ;' -I- 1) and 
per^(M • n) - Hben per^M[fc], with M[b] denoting the submatrix ofM with rows and columns 
labeled by the elements ofb c [n]. 

Theorem 1.2. For any a,peCandA,Be C"^", 

(6) per„(A + B) = J]^ per,(A4 + B40, 

bQln] 

where 4 is a diagonal matrix with {i, i) entry lifieb and otherwise and ¥ is the complement of 
b in [n], and 

n 

(7) per,(AB)= J^ per^(B,) [] A,,,, 

xe[n]" 7=1 

where B^ is the matrix whose jth row is the Xjth row of B and [n]" :- {{ii,...,in) '■ 1 ^ ij ^ 
n for all j = !...,«}. 



SOME ALGEBRAIC IDENTITIES FOR THE a-PERMANENT 3 

Compared to the permanent, the a-permanent has been scarcely studied in the litera- 
ture. At first glance, the a-permanent may not appear as mathematically natural as the 
permanent or determinant: the a-permanent is not an immanant; and it is not clear what, 
if any, interpretation is possible for values of a other than +1. On the other hand, the 
a-permanent arises naturally in statistical modeling of bosons HI [Tsl (called permanental 
processes), as well as in connection to some well-known models in population genetics 
[Tl|. In statistical physics applications, permanental processes are complementary to deter- 
minantal processes, which model fermions in quantum mechanics. The Pauli exclusion 
principle asserts that identical fermions cannot simultaneously occupy the same quantum 
state, which is reflected in the property that detM = if two rows of M are identical. Just as 
the exclusion principle does not apply to bosons, per M need not be zero if M has identical 
rows. 

Because of the practical potential of permanents, devising efficient methods (random 
and deterministic) for approximation is a priority. Recently, there has been some progress 
in this direction ^ |7l |9l; however, a provably accurate method which is also practical for 
large matrices is not yet available. 

In section|2l we prove Theorems ll.ll and ll.2l and observe immediate corollaries. In section 
12.31 we also discuss the relationship betw^een the a-permanent and general immanants. In 
section |3l we discuss computation of permanents in three contexts: in section I3.1[ we 
discuss exact computation of the permanent for some specially structured matrices; in 
section 13.21 we briefly discuss some approximations of the a-permanent based on (O; in 
section I3.3[ we use numerical approximations, see Table [H and inspection of ^ to make 
a conjecture about the computational complexity of the a-permanent, which has not been 
studied. 

2. Identities for the a-PERMANENT 

In addition to (O, (O and ^, we observe several immediate corollaries for the a- 
permanent. We also discuss the relationship between the a-permanent and the immanant 
in section 1231 In section IZTl we prove Theorem ll.lt in section l2!2l we prove Theorem ll.2[ 

2.1. Proof of Theorem ll.li Given n € 'P[n]/ a set partition of [n], we define the nxn matrix 
(uij, 1 < i, i < n), also denoted n, by 

_ J 1, i and / are in the same block of n, 
'^ ■" 1 0, otherwise. 

In this way, the entries of the Hadamard product M • n coincide with the entries of M, 
unless a corresponding entry of n is 0. Since any n, regarded as a matrix, is the image 
under conjugation by a e vS„ of a block diagonal matrix, and the a-permanent is invariant 
under conjugation by a permutation matrix, we can regard M • ti as block diagonal. If we call 
any M • tt a block diagonal projection of M, the Permanent Decomposition Theorem (Theorem 
II. Hi states that the a|S-permanent of a matrix is a linear combination of a-permanents of all 
its block diagonal projections. Also, it should be clear from @ that 



per^(M-7z) = ]^per^M[fc], 



hen 



where the product is over the blocks of n and M{h\ denotes the sub-matrix of M with 
rows and columns indexed by the elements of b c \n\. Because the diagonal product 
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nLi ^j,a(j)'^j,a(f) = for any a whose cycles are not a refinement of n, we have 

n 

per,(M-7T) = 2^a*'^[]M,,(^), 

where, for o £ Sn, we write a < tz to denote that each cycle of a (as a subset of [n]) is a 
subset of some block of n; in other w^ords, a < nii and only if the partition of [n] induced 
by a is a refinement of n. 

To prove (|5]|, we begin w^ith the right-hand side: let a, |S e C and M 6 C"^", then 

fj 

n 

aeSu j=l n>a 
n #a 

aeSu ;=1 ;=1 
n 

aeS„ ;=1 

= Pe^^' 
where s{n, k) := #{partitions of [n] with exactly A: blocks} is the (n, fc) Stirling number of the 
second kind, whose generating function is x" := TJk=i s{n,k)x^'^. Identity ^ follows readily 
by taking a - f^ - -lin ^. 

In the context of point processes, McCuUagh and Moller IITOll show the special case of (O 
with jS = A: € N: 

(8) perteM= Y, ^^*" per«(M • tz), 

w^here 'P[n]:)c is the collection of partitions of [n] with k or fewer blocks. They call this the 
infinite divisibility property of the a-permanent. Also, with M = J„, the nxn matrix of all 
ones, (O relates closely to the Ew^ens sampling formula from population genetics lUIJl. 

More generally, we can take a = -1 in (O and let jS € C be arbitrary, in which case we 
observe 

(9) per^ M = (-1)" ^j (-|S)^*" det(M • tz); 

in particular, for k e IN, 

per_^M= Yj Z:^*" det{M • tt), 

nef[„]:;c 

which plays a significant role in section|3l 
Immediately from ^, we also have 

(-l)"detM= 2]^ (-l/a)^*"per^(M-7r). 
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Previously, Vere-Jones KTSlI obtained an expansion for the determinant as a sum of a- 
permanents, whose form hearkens the MacMahon Master Theorem. 

2.2. Proof of Theorem [Tl 

Proof of identity ^. Given b c [n], we write 4 to denote the nxn diagonal matrix with (i, i) 
entry 



f 1, ieb 
^^^'^' - \ 0, otherwise. 



With this notation, Aly + Bl},c denotes the nxn matrix whose zth row is the /th row of A if 
i Eh or the zth row of B otherwise. 

We prove (O using induction and the cofactor expansion of the a-permanent. For 
M e C"^" and 1 < z, / < n, we write M^''^^ to denote the (zz - 1) x (zz - 1) matrix obtained 
from M by removing its zth row and /th colurrm. For any \ <i <n, the cofactor expansion 
of per^ M along the zth row is 

per^ M = aMi^i per^^ M^'''^ + V Mq per^ M'''^\ 

For n — \,^ clearly holds and we assume it holds for n > 1. 

For convenience, we write M :— A + B and, since it is understood that we are expanding 
along the (zz + l)st row of M, we write M^ := M^"'^^'^\ and likewise for A and B. We also use 
the shorthand [A, B^a,h '■= ^h + BI^, with [A, B]a := [A, Bja^a'^, when fe = fl'^ is the complement 
of fl in [zz + 1]. By induction, we have 

n 

per^ (A + B) = aM„+i,n+i per^M"+^ + Yj Mn+i,j per„ M' 

7=1 
n 

bQ{n\ ;=1 hQMWj] 
n 

foc[n] 7=1 hQMWj] 

n 

foc[n] y=l fcc[„]\{;} 

L perJA,B]b+ J], perJA,B]b 

fcc[n+l]:n+leb fcc[ji+l]:f;+l^b 

D 

Corollary 2.1. For A e C"^" and aeC, 

per„(A + ;„)=2]^a"-*^per^A[b], 

bQln] 
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where J„ is thenxn identity matrix. In particular, 

(10) det{A + I„)=Y^ det A[b] 



bQ[n] 



and 



per (A + I„) = Y\ peiA[b]. 

bQln] 

Remark 2.2. A similar identity to ^ is given in Theorem 1.4 (chapter 2) of IITTll , but, even in the 
case a = 1, that identity differs slightly. We also know, through personal communications, that ([lOl ) 
is known within the context of determinantal processes, and presumably more widely. Regardless, 
we do not formally know of any identities that resemble ^, even in special cases. 

Proof of identity ((7|. To prove (O, we use the cofactor exparision of the a-permanent together 
with induction, as in the proof of (|6]). Let A, B € C"^" and a € C. For n = 1, (O is trivially 
true. In the induction step, we assume (O holds for n > 1. 

For X = (xi, . . . , Xn) e [n]", we write Bx to denote the n x n matrix whose (/, /) entry is 
Bx(/, /') := B(x„ ;'), i.e. the rth row of Bx is the x,th row of B. We write M := AB and, again, 
M^ := M^"'^^'^^ denotes the submatrix of M obtained by removing the n + 1st row and /th 
column. We have the following for A,B e c("+i)><("+i). 

n 

per^(AB) = aMn+i,n+i per^ M"+^ + ^^ Mn+ij per^ M' 

n n n 

xe[n]" k=l ;=1 xe\n\" 



Xt 



k=l 



fn+1 



n fn+1 



k=l 



, i=l J xe[n]» k=l j=l \ i=l J xeln]" k=l 

( n+1 \ n 11+1 n n 

Y^ A„+i,B,,„+i Yj p^'a sf n ^'^'-'c + L ^"+1'' L ^i'j L p^'a sr' n ^^'^ 

) xe[n]" k=l i=l j=l xe[n]" 
n+1 n 

/=1 ;=1 

H+l n+1 

= Y L Per«B.n^^'^ 



n+l 



xe[n]" 



)c=l 



=1 xe[n+l]"+i:x„+i=i 



*:=1 



n+l 



= Y P^'a^^Il^k,x,. 

.Te[n+l]«+i ^=1 

This completes the proof. n 

Remark 2.3. Theorem 1.3 (chapter 2) in Mine [TT\ restates a theorem ofBinet and Cauchy for the 
permanent of a product of matrices, but that identity is different than (|7|. Both ^ and (O are more 
suitable to statistical applications than their counterparts in [llj, because they give the explicit 
normalizing constant of a probability distribution on subsets of[n], which is relevant to various 
applications involving clustering. 
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2.3. Immanants and the a-permanent. In representation theory, the permanent corre- 
sponds to the immanant associated to the trivial representation of Sn, while the deter- 
minant is the immanant corresponding to the alternating representation, the only other 
one-dimensional representation of Sn- With A -\ n indicating that A is a partition of the 
integer n, the immanant indexed by any A H n is 

n 

(11) ImA M-.^Yj '^'^(^) n ^^■''^0> 

aeSi, 7=1 

where XAio) is the character associated to the irreducible representation S^, the Specht 
module corresponding A -\ n; see [3]. Aside from a - +1, the a-permanent @ does not 
correspond to an immanant of a particular index. However, there is a precise relationship 
between (HJ and ((TT]l . which we now discuss. 

2.3.1. Characters and immanants. For each Ann, the character xa : 5„ ^ R associated to 
the irreducible representation p,\ is a class function, i.e. XAio) = xa{t^ot~^) for all o,t e S„. 
The conjugacy classes of Sn correspond \.of„:-{A-\n], the collection of integer partitions 
of n, 

Tn ■■= {{ni, . . . ,nk) : ni > ■ ■ ■ > nk > 0, ni + ■ ■ ■ + nk ^ n] . 
Moreover, the collection (xa, A 6 'Pn) is orthnormal and forms a basis for the space of class 
functions on S„. The a-permanent is a polynomial in a and, though a** is not a character, 
it is a class function. Hence, for every a there exist constants {c", A e Pn) such that 

(12) per^M = 2]^c«ImAM. 

A-\n 

Special cases are a = +1, for which c^^ = 6A,(fj)/ the Dirac mass at A = (n), andc^^ = (-1)"6a,1"/ 
the Dirac mass at A = 1" := (1, . . . , 1) e Pn- 

For n 6 IM, let X := (Xa(v), A, v h n) be a Pn ^'Pn matrix, i.e. its rows and columns are 
labeled by Pn, with (A, v) entry equal the value the A-character takes on the conjugacy class 
indexed by v. For any (A, v) 6 !P„ x Pn, we write X^'^'^'^ to denote the minor determinant of 
X with A row and v column removed. We then define Y := (X''^'^', A,v -\n) as the matrix of 
minor determinants. 

Theorem 2.4. For every A e Pn, define the function ca : C — > C fey 

(13) CA(a) := ^ y a*"XA(cT). 

n! '—^ 

aeS„ 

Then the a-permanent satisfies 

per^M= Yj CA(a)ImAM 
AePn 

for all a eC and M e C"^". Alternatively, c := (cA(a), A -\n) solves the equation 

(14) A = Xc, 

where X is defined before the theorem, A := (a*'^, A -\ n) is a P„ xl column vector and #A denotes 
the number of parts of A. The solution to ([14]) is 

(15) c = YA/detX, 
where Y is defined above. 
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Proof. That such a ftinction c;^ exists is a consequence of ((T2|) and the surrounding discussion. 
To obtain l(T3]l , we note that 

n n 

per„ M = 2]^ a*" [] M^,,(y) ^ YjTj ^'^(«)^^('^) R ^/''^O) ^^^ ^^ ^ ^ C"^"; 

whence, 

(16) a*'' = Yj Ca (a);rA (c^), for all oeSn- 

A-\n 

The Fourier inversion theorem gives ((T3)) . 

To see ((15)) , it is enough to show that YX = det(X)7f>„, where 7f>„ is the Pn x "Pn identity 
matrix. Indeed, by orthogonality of characters, 

(YXW = E..(.)X<-> = {"-". J^^^ 

for all A, V H n. Orthogonality of characters, and hence the columns of X, implies detX i^ 0, 
completing the proof. D 

Corollary 2.5. For any a,p eC, 

(17) per^ M = Y, c/ Y. (-a/|3)-L#" ImA(M • tz), 

where c^ := Ca(-j6) zs defined in ([l3l ). 

Proof. Immediately, by combining (|5ll and ((12)) , we obtain the a-permanent as an explicit 
linear combination of immanants, in terms of the coefficients (ca(jS), A H n), for arbitrary 
)3 € C, yielding ([17]). □ 

From ((TT)) , we can write 

n 

per^(M • 71) = 2]^ n Mj,a(j) Y^ c{xa{o) 

a<n j=l A^n 

For a e Sn and n 6 ^[„], we write a ~ n to denote that the cycles of a correspond to the 
blocks of 7i; and, since Xa is a class function, we define XAi^i) to be the common value that 
Xa takes on [a ~ ti}. Defining f{n) = per^(M • ti) and g{n) = L„^„ n^li Mj,a(j) T^x^n <^xXA{n), 
we obtain 



(18) 



n 

a~n j=l 



vHn (7e5„ Tc'<n ben 



as a consequence of the Mobius inversion formula ( 1(161 , proposition 3.7.1 and example 
3.10.4). As a special case, when M = /„ in 



[](#& - 1)! L L /^*'';rv(TT);rv(a) = Y n /^'*' n(-i)'"''"'(#"[^ - !)'• 

foSTi vHii 0g5„ n'<nben' ben 
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3. Computing the a-PERMANENT 

3.1. Permanents of special matrices. Although the a-permanent is difficult to compute in 
general, it can be computed explicitly for some specially structured matrices. For example, 
for any a € Sn, let Pa be its matrix representation. Then per^P,7 = a*'^. Also, if, as before, 
we regard n e !P[„] as a matrix, then 



per^Ti 



n 

ben 



a' 



the product over the blocks of the rising factorial a^*^ := a{a + !)■ ■ ■{a + #b-l). Finally, for 
a 2 X 2 matrix 

...Mil ^12\ 
\A21 A22I' 

define A[ni, n2] as the (wi + ni) x (nj + n2) block matrix with (i, j) entry 

' All, i,i<ni 

AV !/• -S A22, i, i > «1 



A 



21/ 



]>ni 
i > ni, i < ni. 



For example. 



A[2,3]- 



'All All A12 A12 Ai2^ 

All All A12 A12 A12 

A21 A21 A22 A22 A22 . 

A21 A21 A22 A22 A22 

.A21 A21 A22 A22 A22, 
Such block matrices arise naturally in statistical modeling, and the follow^ing identity w^as 
obtained in [14J. 

Proposition 3.1 (Proposition 1 in A.2 [14]). Let A be 2x2 as above and define 

A11A22 



A12A21 



Then 



"4^^^ n%^^pi 



jlaV 



per^A[ni,n2]=A';|A^>Tni^Tn2 J]^ J. 

7=0 

where a^"^a^"^/a^> := whenever the quantity is undefined. 

Another natural consideration are what we call homogeneously symmetric matrices H[fl, b; n\. 
For a,b eC, we define H[fl, b; n] to be the nxn matrix with diagonal entries a and off-diagonal 
entries b. Since n is typically understood, we will usually write H[a, b]. For example, 

fa b b\ 



H[a,b;3] 



b a b 
b b a 



Proposition 3.2. Let H[fl, b] be a homogeneously symmetric matrix and write d := b/a. Then, for 
any a eC, 



(19) 



per„ H[fl, fo] = «" y y cin, k, l)a''d"-\ 



k=l /=0 
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where c{n,k, I) denotes the number of permutations ofn having exactly k cycles and exactly I fixed 
■points. 

Proof. It is plain that H[fl, &] = aH[l,b/a] and per^^H[fl, fo] = fl"per^H[l,d], d :- b/a. So it 
is enough to consider matrices with unit diagonal. Identity ( fT9l ) follows directly from the 
definition of the a-permanent and the coefficients c{n, k,l). n 

The numbers c{n,k,l) are related to the rencontres numbers f{n,l), where f{n,l) is the 
number of permutations of [n] with exactly / fixed points. The rencontres numbers corre- 
spond to the number of partial derangements of an n-set, with a derangement defined as a 
permutation with zero fixed points. By this definition. 



f{nJ)^Yj'^n,k,ll 

k 

prompting us to call c{n, k, I) the (n, k, l)-generalized rencontres number. The rencontres num- 
bers satisfy /(0, 0) = 1, /(1, 0) = and, for n > 1, 

f{n + l,0) = n{f{n,0)+f{n-l,0)). 

For <l <n, f{n, I) satisfies 

f{nj) = i^^fin-l,0). 

The (n, k, I) Rencontres numbers satisfy 
(20) 

k i)-i cin - l,k - 1,1 - 1) + {n - I - l)c{n - l,k,l) + {I + l)c{n - l,k,l + 1), <l <k <n, k <n-l 
1 0, otherwise. 

The recursion (|20)) has a simple combinatorial interpretation. We call a € Sn an {n,k,l) 
permutation if it has exactly k cycles and exactly I fixed points. An {n, k, I) permutation can 
be obtained from o € Sn-i by inserting n + 1 into one of the cycles of o only if 

(1) a is an (n - 1, / - 1, A; - 1) permutation: in this case, we obtain an (n, k, I) permutation 
by appending n + 1 to a as a fixed point; 

(2) a is an (n - 1, k, I) permutation: in this case, we insert n + 1 into one of the k — I cycles 
that are not fixed; there are n — I — 1 possible points of insertion; 

(3) o is an {n — l,k,l + 1) permutation: in this case, we choose from one of the / + 1 fixed 
points and insert n + 1, bringing the number of fixed points in the new permutation 
to I while keeping the total number of cycles fixed. 

The recursion ( [20l ) makes these numbers easy to compute. The {n,k, /)-generalized rencon- 
tres numbers are also related to the number of derangements. Let g{n, k) denote the number 
of permutations of [n] with k cycles and zero fixed points. Clearly, gin,k) = c{n,k,0) and 
/(n,0) = T,kS(^'^)- ^6 ^Iso have the expression 

where 5ki - 1 H k - I and otherwise. Appendix 14.21 contains tables of the generalized 
rencontres numbers for n < 10. 
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3.2. Approximating the a-permanent. This and the following section concern the approx- 
imation and computational complexity of the a-permanent. Among previous algorithms 
and asymptotic approximations of the permanent, the randomized method in |6| seems the 
most practical, but is not provably good. Other methods, such as those in ||5l[Zl, are either 
not practical (the method in |j5l still requires 0(n''(log n)^) operations) or not accurate (the 
approximation in [7] is accurate up to a factor e"). All of the aforementioned approximation 
methods are valid for non-negative matrices and, in the case of |6|, a > 0. The assumption 
that all entries of M are non-negative is, in some sense, the permanental analog to positive 
semi-definiteness in the determinantal case. For a > and M,y > 0, each term in H]) is 
non-negative. Analogously, for M positive semi-definite, det(M ■ n) >0 for all n e P[„] by 
Sylvester's criterion. 

Positivity of the terms of @ is crucial to the approximation in |6|, which uses the 
statistical method of importance sampling (IS) to approximate per^ M. An introduction to 
IS is given in [13J pp. 411-412. In IS, we estimate the sum 

n 

aeS„ i=l 

by drawing random permutations ai, . . . , aw, for N as large as is computationally feasible, 
and, writing f{a) := a**^ DLi ^j,a(j), we approximate per^ M by 

w^here P{a) is the probability of a. The choice of per^M is so that our estimate is unbiased and, 
by the law of large numbers, if N is large enough, per ^M should be close to the true value 
per^ M. The efficiency of this method relies on choosing P which is close to optimal. It is 
known that P{a) oc f(a) is the optimal choice, which is rarely practical. Kou and McCuUagh 
use a variation of importance sampling, called sequential importance sampling, which 
allows them to draw from P which is nearly optimal. 

Because the determinant can be computed efficiently, (|9]l is a natural candidate for use 
in IS. To use (|9]) in IS, we choose a sequence 7ii, . . . , tin of random partitions from P* and 
estimate (-1)" per^ M by 

( 1^"— AA V (-«)^*"det(M-7i) 
(-1) per^M = l^ — . 

7=1 ^ 

In view of Q, the optimal importance sampling distribution is P*{n) oc K-a)-!-*" det(M • n)\. 
One issue is that, even if we could draw from P* and we assume M is positive definite, for 
any a e IR, there will be both positive and negative terms in Q, and our estimate based on 
importance sampling could have extraordinarily high variance. In fact, some preliminary 
numerical tests illustrate this issue, but suggest that when a is a negative integer, the 
problem is simplified and a reasonable approximation is possible. 

3.2.1. Some numerical illustrations. This section is a prelude to section [331 where we draw on 
observations from numerical tests to analyze, intuitively, any further reaching implications 
for the complexity of the a-permanent. We show only a few numerical approximations in 
order to illustrate IS, as well as some of the caveats discussed above. 
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Actual 


Estimate + standard error 


relative error 


per_2 Xi 


407.52 


406.43 + 10.32 


2.53% 


per_3 Xi 


117488 


117300 + 3957 


3.37% 


per_2 5 Xi 


-44088 


-427609 + 7.03 x 10^ 


1643% 


perXi 


1.6 X 10^ 


3.1 X 10^ + 4.8 X 10'^ 


1548% 


-2-permanent 




199.19 + 8.10 


4.01% 


-3-permanent 




41692 + 1876 


4.45% 


-2.5-perinanent 




43338 + 995181 


2296% 


1-permanent 




-1.1 X 10^ + 6.3 X 10^ 


5727% 



Table 1. Estimates of the a-permanent for a = 1, -2, -2.5, -3 for randomly 
generated, symmetric positive definite matrices. The first four rows are the 
a-permanent for Xi, given in the appendix. Each of the bottom four rows 
reflects 200 iterations of generating random, symmetric positive definite 
matrices (8 X 8) and comparing, for each iteration, the estimated quantity 
and the actual quantity. 



The sampling distribution P* we used is the Pitman-Ewens(fl, 0) distribution [12], for 
various choices of (a, 0). We would not expect that this distribution is near optimal; however, 
when a — —k for k e IN, we can choose a < and - -ka, which restricts the Pitman-Ewens 
distribution to P[n]:k- On the other hand, when a is not a negative integer, it seems there 
is no optimal choice of {a, 6); in this case, all estimates, based on different choices of {a, 6), 
returned values similar to those in Tabled] for the 1 and -2.5 case (the table reflects a — 
and = 1 in these cases, the standard Ew^ens sampling formula [2J). The purpose of the 
demonstration is to show the striking difference when using (|9]) and a naive choice of P* to 
approximate the a-permanent in the cases when a is a negative integer and all other cases. 
The summaries are given in Tabled] 

The standard errors for estimating the ordinary (a = 1) and -2.5-permanents exceed 
the value of the quantity being estimated. In these cases, our approximation is unstable 
because it involves a sum of both large positive and large negative numbers, a consequence 
of the coefficient (-jS)-'-*" in (O, which is alternating between positive and negative values 
for #n > a. Stability is only achieved for negative integer values. 

3.3. Computational complexity of the a-permanent. We conclude with some brief re- 
marks about computational complexity of the a-permanent, which, to our knowledge, has 
not been studied. For any a € R, we have 

per^M = {-!)" ^j (-«)^*" det(M • ti). 

It seems intuitive that per^ M for a > can be no easier to compute than per M; however, 
when kis a negative integer, we have cancellation of all terms indexed by n with #n > —k, 
yielding 

per_jtM= Yj /c^*" det(M • 7z). 



neP, 



l,l]:k 



When M is positive semi-definite, all terms on the right-hand side above are non-negative 
and it seems reasonable that computation might be easier than the general a case. 
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Even if M is not positive semi-definite, a reasonable approximation might be possible 
because the sum is over 'Piny.k, which is asymptotically of much smaller order than P[n]- 
Writing B{n,k) := {tt € ^[„] : #n = k}, the {n,k)th Bell number, and B{n) :- Lit^C"/^)/ we see 
that, for fixed fc 6 N, 

k 

B{n,<k):=YjBin,k)-'ie'/k\ asn^oo 

7=1 

and B{n) ~ „-i/2[^(„)]«+i/2gA(n)-,i-i as n ^ oo, where A(n) = e^^"\ for W{n) the Lambert 
W function. Identity (HJ is also consistent with our knowledge that computation of the 
determinant is much easier than the permanent; the case j6 = -1 is the degenerate sum over 
one term in ^. Our investigation prompts the follow^ing conjecture. 

Conjecture 3.3. Computation of the a-permanent is #P-complete, except for a e {-1, ... , -K},for 
some _K £ N U {oo}. The complexity for a € {-2, -3, . . .} lies somewhere between P and #P and is 
monotonically non-decreasing as a decreases. 

We leave open the possibility that X = 1 above, in which case our conjecture is that the 
a-permanent is #P-complete except for a - —1. The crux of our conjecture, however, is that 
the a-permanent must be (at least) #P-complete unless a is a negative integer. 

4. Appendix 

4.1. Matrices for numerical approximation. Below is the matrix Xj used in approxima- 
tions in TablelH The matrix is symmetric about the diagonal and we include only the upper 
triangular part of Xj . 

^4.42 



Xi 



3.13 3.14 


3.45 


4.01 


3.85 


3.39 


2.70> 


2.70 1.99 


2.44 


3.07 


2.83 


2.27 


1.84 


3.52 


2.26 


2.73 


2.43 


2.85 


2.36 




3.57 


3.01 


3.17 


2.93 


1.90 






4.12 


3.69 


3.01 


2.03 








3.91 


3.03 

3.27 


2.08 
2.22 
2.33, 



4.2. Generalized rencontres numbers. Here, we list tables of the generalized rencontres 
numbers from section 13^1] for n < 10. Blank entries correspond to 0. 



(2X1) 


1=0 1 2 


k=l 

2 


1 

2 1 



(3,fc,/) 


12 3 


1 

2 
3 


2 

3 

1 



(5,^,/) 





1 


2 


3 4 5 


1 


24 








2 


20 


30 






3 





15 


20 




4 











10 


5 











1 



(4,^,/) 





1 


2 


3 4 


1 


6 








2 


3 


8 






3 








6 




4 











1 



(6,fc,/) 





1 


2 


3 


4 5 6 


1 


120 










2 


130 


144 








3 


15 


120 


90 






4 








45 


40 




5 














15 


6 














1 
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(7,fc,/) 





1 


2 


3 


4 


5 6 7 


1 


720 












2 


924 


840 










3 


210 


910 


504 








4 





105 


420 


210 






5 











105 


70 




6 

















21 


7 

















1 



{8,Kl) 





1 


2 


3 


4 


5 


6 7 8 


1 


5040 














2 


7308 


5760 












3 


2380 


7392 


3360 










4 


105 


1680 


3640 


1344 








5 








420 


1120 


420 






6 














210 


112 




7 




















28 


8 




















1 



(9X1) 





1 


2 


3 


4 


5 


6 


7 8 9 


1 


40320 
















2 


64224 


45360 














3 


26432 


65772 


25920 












4 


2520 


21420 


33264 


10080 










5 





945 


7560 


10920 


3024 








6 











1260 


2520 


756 






7 

















378 


168 




8 























36 


9 























1 



{W,k,l) 





1 


2 


3 


4 


5 


6 


7 


8 9 10 


1 


362880 


















2 


623376 


403200 
















3 


303660 


642240 


226800 














4 


44100 


264320 


328860 


86400 












5 


945 


25200 


107100 


110880 


25200 










6 








4725 


25200 


27300 


6048 








7 














3150 


5040 


1260 






8 




















630 


240 




9 


























45 


10 
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